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We consider the simplified dielectric breakage model used for simulation of the kinetics of cellular 
material breakage under the pulsed electric field (PEF) treatment. The model is based on an 
effective media approximation, which includes equations with the same morphology parameters 
as in percolation theory. The probability of a whole cell breakage by the pulse with ti duration 
is estimated on the basis of electroporation theory. We account for the bridging effect resulting 
1-^ ' from the deviations of the local conductivity near the selected cell from the average effective media 

conductivity. The most important feature of the proposed model is the existence of the "jamming" 
behaviour occurring sometimes in experimental observations of the biological tissue breakage. The 
different transitions corresponding to the "jamming" steps are identified. The experimental results 
are obtained for thin apple slices treated with electric pulses at field strengths E — 0.2 — 2.2 kV 
2 ' cm~^, pulse durations ti — W — 100 fis, pulse repetition times t = 10 — 100 ms and the number 

C/3 ' of pulses A*' = 1 — 100000. The model gives results consistent in general with the experimental 

observations. We discuss the correlation between the degree of cellular material destruction, field 
strength, time of PEF treatment and power consumption. 
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Q ' Notation 

Cm specific capacity of membrane, F m^- 

^ ^ C* = C„(e^/e™-l)/(27) 

►^ I C specific heat of cellular material, J kg^^ K"-*^ 

dm membrane thickness, m 

00 ! c^c cell diameter, m 

E electric field strength, kV cm"-'^ 

^«0 ' AF* = t:lu^ / (kTj) , reduced critical free energy of pore formation 

g breakage probability for a membrane 

j current density, A m.~^ 

k Boltzmann constant, 1.38065812xl0"23j K^^ 

j^ ■ 1/ = dcNz, total width of a sample 

m = dc/dm 

n number of intact cells neighbouring of a given cell 

'^ ' iV number of pulses 

G , Nb number of pulses at the moment of the total dielectric breakdown of material 

O ' Nm number of intact cells 

Nx,Ny,Nz dimensions of a lattice 

P degree of biological tissue destruction 

P„ local degree of cell destruction 

. . . Pp critical percolation degree of destruction 

C^ ; rm,c.i = dm,c,c/iSm(Tm,c,i), resistance, Ohm 

R total resistance of a linear chain. Ohm 

Sm '^ d^ , cross section area of a cell or a membrane, m^ 
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t scaling exponent 

ti pulse duration, fis 

At pulse repetition time, ms 

T temperature, K 

AT temperature increase over initial value, °C 

u transmembrane voltage, V 

Uo midpoint of a probability transition function g, V 

Au width of a probability transition function g, V 

U external voltage, V 

W moisture content, % 

Greek letters 

a geometry factor 

7 surface tension of membrane, N m"-'^ 

Ew =80, dielectric constant of water 

Em = 2, dielectric constant of membrane 

^ = P/(l - P) 

A = dcdrn/dmO'c 

p density of cellular material, kg m^^ 

(7 conductivity, Sm~^ 

(Ti = (7c\/{l + X)^ conductivity of an intact cell 

2-1 = Zje/zjc 

V — 1/* 

^e,c,i — ^e,c,i 

T lifetime of a membrane 

Too parameter, lifetime of a membrane at T = oo 

V' = i/Pp - 1 

Lo linear tension of membrane, N 

Subscripts 

a after treatment 

b before treatment 

c cell 

e effective 

i intact cell 

j juice 

m membrane 

t total 

Abbreviations 

PEF pulsed electric field 

LC linear chain 

MF mean field 

RSA random sequential adsorption 

RSO random sequential occupation 

I. INTRODUCTION 

Pulsed electric field (PEF) treatment is a new and promising nonthermal processing method for heterogeneous 
cellular materials. PEF methods are based on the effect of cell transformation or rupture under an external electric 
field, which results in increase of the electric conductivity and permeability of cellular material. This effect known 
as dielectric breakdown (Zimmermann, Pilwat & Riemann, 1974), or electroplasmolysis (Scheglov, 1983), can be ex- 
plained generally by electroporation, i.e., electric field induced formation and growth of pores in biological membranes 
resulting from their polarization. The method of electroporation became very popular because it was found to be 
an exceptionally practical way of transferring drugs, genetic material (e.g. DNA), or other molecules inside the cells 
(Chang, Chassy, Saunders & Sowers, 1992). Recently, scientists began to use the PEF methods for treatment of liq- 
uid food (fruit juices, milk etc.), as an alternative to high temperature preservation (Barbosa-Canovas, Pothakamury, 



Palou & Swanson, 1998; Hulsheger, Potel & Niemann, 1983). PEF protocols also have been recently introduced for 
cellular tissue treatment (Knorr, Geulen, Grahl & Sitzmann, 1994; Knorr, 1997). 

However, presently the wide spreading of PEF treatment for the purposes of nonthermal processing of heterogeneous 
cellular materials is restrained by the poor reproducibility of experimental data, the unclear mechanism of dielectric 
breakdown in the cellular systems, and the absence of criteria for choosing optimal parameters of HEFP treatment. 
The most optimal for the tasks of this kind is to apply computer simulation methods. There exist many types 
of models for simulation of dynamic breakage processes in random disordered systems. The most popular is the 
model of random resistor (or random conductor) networks (Arcangelis, Redner & Herrmann, 1985; Duxbury & Beale, 
1995). The models of this type are applied for simulation of the physical behavior of granular superconductor, metal- 
insulator composites and various other disordered materials. Usually, the breakdown process kinetics is investigated 
at the network models of different types through the iterative solving of the Ohm-Kirchhoff 's system of equations 
with allowance for boundary conditions and with choosing the relevant initial conditions for breakdown initialization 
(Lebovka & Mank, 1992). But all the above mentioned models are highly time consuming and they require a lot of 
CPU-time for the real cellular material simulation (Press, Teukolsky, Vetterling & Flannery, 1997). 

The objective of this study is to develop the simulation model describing the breakage kinetics under conditions 
of PEF treatment and to correlate experimental data for heterogeneous cellular materials with simulation results. 
We consider the very simplified model of cellular material breakage under the PEF, which is based on a generalised 
approximation of effective medium. This model takes into account the percolation properties of medium and local 
bridging effects. This approach is, certainly, rather rough, same as all the other effective medium approaches, but it 
allows us to obtain the realistic description of breakage processes and trends of kinetics. 

II. BACKGROUND THEORY AND SIMULATION MODEL 

A. Field Induced on a Cell Membrane in an External Field 

We consider at first the case of an individual cell inside the homogeneous medium having a low effective conductivity 
ae placed into an external electric field of strength E. The maximum voltage u between the internal and external 
cellular surfaces of membrane (or transmembrane voltage) induced on cell poles by the external field, is equal to 
(Neumann, Sprafke, & Boldt, 1992): 

u = afdaE (1) 

where a is the geometry factor equal to 0.75 for the spherical geometry of cell and to 1 for a cell of cubic geometry, 
and / depends on the electrical and geometrical properties of the cell (Kotnik, Miklavcic & Slivnik, 1998) 

„ _ 0-e(6TO0-e - (12t?7? - 8m^)(fT^ - (Te)) , , 

{a„i + 2ae){am + <Jc/2) - (1 - 2m)^{ae - (Tm){iyc - (Jm) 

Taking into account the low conductivity of membrane, Uml Oc <<1 and expanding Eq. (0) in power of crm/crc, we 
obtain: 
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I = {l + ^)/K + 0[{a„Ja,Y\, (3) 

where K ^ l-\- A(2 + (Tc/ce)/4 , and under the usual conditions (standard values of dc, dm, Cm, Cc and ae parameters 
for this problem are collected by Kotnik, Miklavcic and Slivnik (1998)) A ~ 10^^ — 10^^ << 1 and / ^ 1. For 
the purposes of estimation we shall consider further on only the cells of cubic geometry, so we put a ~ 1. In this 
approximation, / is equal to the normalised transmembrane voltage u/ d^E. 

Now we consider a problem of transmembrane voltage u estimation for a cell inside a biological tissue. In this case, 
the external media is not homogeneous, so, Eqs. (El) and (0) are inaccurate. They can be used only as approximate 
estimation by assuming a certain effective value of Ge- This is just the mean field (MF) type of approximation, which 
does not take local conductivity fluctuations into account. It can be assumed, as a rough approximation, that in an 
inhomogeneous medium the equation (1^) may also be used for determination of the u values, however, / = /(r) will 
depend on the space co-ordinate r and will vary from cell to cell. 

As it follows from Eq. (H), the value of f gradually increases with the increase of external medium effective 
conductivity a^, which depends on the biological tissue destruction degree, P. Here we define the biological tissue 
destruction degree as: 



P=l- N,n/Nt. (4) 

As a more precise estimation of cr^ versus P dependence, we have used the generahsed effective media equation 
with the same morphology parameters as in percolation theory (McLachlan, 1989): 

1 + V'E +E„ + ^I]-"- ^^> 

This equation reduces to the well-known Bruggeman's symmetric equation when i = 1 (Landauer, 1978). Every- 
where throughout this work we use the values oi t = 2.3, Pp — 0.307, which are typical for the 3d simple cubic lattice 
(Sahimi, 1993). In the general case we obtain the following solution of Eq. (g): 



<7e{P) = <yc 



ipjijj - T,o) + Y.qiIj - 1+ \/(9!( ,i\\ 

+ ^ipH^P + So)2 + 2(^[Eo(^ + 1)2 - 7A(So - 1)2] + (^Eo + 1)2 ' / ^^^^'^ + '>> 



(6) 



which can be easily used for the numerical estimation of a^ versus P dependency. 

From Eq. (g) we can easily obtain a^ = fJi = ac\/{l -I- A) at P = and ae—Uc at P = 1. 

The example of / versus P dependence on the mean field approximation obtained with the help of Eqs. (0) and (o) 
is presented in Fig. |l| (line 1). At low degrees of the biological tissue destruction, P ^ 0, the effective conductivity 
approximately equals to Ce = Ci = (TcA/(1 -I- A). At complete destruction of a biological tissue, P ^ 1, we have 
(7e ~ <7c. Thus, we obtain the following limiting relations in MF approximation: 

/(P^0) = 1/(2 + 3A) « 1/2, and (7) 

/(P^1) = 1/(1 + 3A)«1 

i.e., the value of / approximately doubles with the biological tissue destruction increase within the range of P = — 1. 
In fact, the change in / with increasing biological tissue destruction degree (P) is impossible to describe within 
the frames of such simplified approach, because of substantial local fluctuations of the external medium conductivity 
existing in the real inhomogeneous biological tissues. Such fluctuations are especially high near the punctured cells, 
where the local values of conductivity can exceed the effective average value a^ by some orders of magnitude. The 
MF approximation works especially poor at low degrees of biological tissue destruction, P — > 0, when the local 
conductivity fluctuations are high. In order to demonstrate the complex behaviour of / versus P dependence in a 
real biological tissue, we have considered the oversimplified linear chains (LC) approximation as a simulation of a real 
tissue structure. The tissue is simulated by linear chains consisting of Nt cells of a cubic-like geometry (a = 1). The 
total resistance of such linear chain consisting of serially connected membrane resistances (r^ = dm/icmSm)) and 
intercellular fiuid resistances {vc = dd {(JcSm)) can be estimated using the following relation: 

R^{Ntr,+Nrarra)- (8) 

The density of current running through such linear chain is equal to: 

J - U/{RS^) = U/{NtdJa, + Nrad^/a^) = a^U/{Ntd,) = a^E, (9) 

where 

ae = Aa,/(A+l-P), (10) 

is an effective chain conductivity, and E = U/L = U/{Ntdc). 

The percolation point is observed in LC-approximation at P = Pp = 1, and it can be easily shown, that the more 
general Eq. (||) reduces to the Eq. ( [10| ) in the case of t = 1. 

The voltage drop at a single membrane is equal to: 

U = jSmTm = jdm/cTm = dcE(7e / {Xac) , (11) 

or we obtain for / = u/dcE (at a = 1): 

f^a,{P)/\a,. (12) 

The example of / versus P dependence for the linear chain approximation obtained with the help of Eqs. (0) and 
(H) is presented in the Fig. |l| (line 2). 



We obtain from Eqs. (|6|) and ( ]12| ) in LC approximation with account for A << 1: 

/(P^O) = 1/(1 + A) «1, and (13) 

/(P^l) = l/A»l 

i.e. in this case the value of f can considerably increase in the course of biological tissue destruction, and it is in 
contradiction with the conclusion, which follows from Eq. (|^) on the basis of MF approximation. 

Moreover, MF approximation gives us f{P — > 1) « 1, but in LF approximation we have f{P ^ 0) = 1/(1 + A) w 1. 
How we can clear this contradiction? MF approximation is good only at large values of P — > 1, when all the cells 
of tissue are destroyed and the media is practically homogeneous. On the other side, LC approximation works well 
only at small values of P — > 0. In this latter case we can neglect any bridging effects, which are very important in 
the middle of the range of P values. The bridging effect results from the deviations of local conductivity near the 
selected cell from the average effective conductivity of the whole media. Analysing the parallel electrical circuits in 
the near-neighbours environments of a given cell, we can easily approximate the averaged local cell resistance r„ as 

r„«l/((l-P„)/r, + P„/r,). (14) 

For the cells arranged in the sites of a simple cubic lattice we have P„ = (1 — n/6), where n is the discrete number 
of intact cells neighbouring a given cell. Then the voltage drop at a single membrane equals to: 

u = jrnSmrm|r^ = - — \ ^ , (15) 

and we obtain the following generalised equation for estimation of the local value /„ = u/dcE 

/„(P, P„) = aAP)/{X + PnW (16) 

where in the general case the value of ae at given P may be determined from Eq. (|6|). 

Figure shows the example of local conductivity /„ versus P dependencies obtained at various n with the help of 
Eqs. (H) and (|lq). We see that the local behaviour of fn{P,Pn) function can be rather complex and it can reflect 
the local conductivity changes near the chosen cell. We have shown with the dashed curve the imaginary fluctuations 
of /„ with P increase and successive breakdown of neighbouring cells. We understand that this approach proposed 
for estimation of the /n(P, Pn) behaviour is very simplified. Particularly, here we overestimate the importance of the 
lattice contribution to the /„ (P, P„ ) behaviour due to the discreteness of this problem, and in the case of continuum 
more smooth behaviour of fn(P, Pn) can be observed for the continual parameter P„. 

The behaviour of the averaged value /(P) =< fn(P,Pn) > (here <. . . > implies averaging over all the cells of a 
system) will depend on the spatial distribution of broken cells in such system. Let us consider as an example, the 
/(P) behaviour for the model of random sequential occupation (RSO) of the lattice sites with the broken cells. This 
problem is very complex and has no exact solution. Here, the Monte Carlo simulation is a rather simple and useful 
method allowing to get the /(P) function for the given type of cell distribution in the system. We have considered 
the cells located in the sites of a simple cubic lattice with the following dimensions: 100x100x100. The simulation 
consisted of successive random choice of an intact cell, its further destruction and averaging of the calculated /-values 
over the whole lattice. We used the periodical boundary conditions along all the x, y and z directions in order to 
reduce the influence of boundaries. The results were averaged over 10 different initial configurations. The example of 
/ versus P dependence obtained as described above is presented in the Fig. ^ (line 3). The dashed line 4 corresponds 
to / dispersion (A/) and characterises fluctuations of / in the system. We see that on increase of the destruction 
degree (P) the values of / first decrease, then pass through a minimum and increase again. The point of minimum 
/(P) approximately corresponds to the maximal fluctuations of the local values /„ in a system, i.e., to the maximum 
of A/. The important positive feature of this method of /(P) estimation is as follows: we get / w 1 in both limit 
cases of P — > and P ^- 1, i.e., we remove here the contradiction between MF and LC approximations. 

B. Probability of a Single Cell Destruction 

The average lifetime of a membrane in the external electric field can be estimated with the help of the following 
equation (Weaver & Chismadzhev, 1996): 

r(T,u) = rooexp(AP7(l + M2C*)). (17) 



Lebedeva (1987) has presented the following estimates for the general lipid membranes: Tqo ~ 0.37x10 ^ s , 
Lo ^ 1.69x10"" N, 7 « 2x10-3 N m'^, Cm ~ 3.5x10^^ F m'^ at T = 298 K. 

Then the breakage probability for a membrane (as a whole cell) during the impact period of an impulse with 
duration of ti may be estimated as 

g = l~exp(-i,/T(r,u)). (18) 

Taking Eq. (|l7| ) into account, we can rewrite Eq. ( |l8| ) in the following convenient dimensionless form 

g{u') = l-exp(-ln2/expa([l - (1 - M*2)/(aAuln2/uo)]"^ - 1)), (19) 



where u* = u/uo, Uo = y^{AF*/a - 1)/C*, Au = Uo ((1 - a/AF*)aln2)" , a = ln(ti/(Too ln2)), and du > duc = 
l/(aln2) . 

We see that g{u*) is a kind of probability transition function and u = Uq corresponds to the midpoint, where 
g{u) = 1/2 and du is the width of this function. We can estimate the following actual values for the general lipid 
membranes using the experimental data of Lebedeva (1987): Uo ^ 0.92V, Au ~ 0.41, a ~ 3.66 and Auc ~ 0.394 (at 
t, = 10 /is) and Uo ~ 0.71V, Am « 0.26, a « 5.97 and Au^ ~ 0.242 (at t, = 100 /xs). 

C. Description of the Simulation Model 

To achieve the best performance, we have used here a hybrid of MF and cellular automaton strategies in order 
to reduce the computational complexity. First we construct an array of undestroyed cells located in the sites of a 
simple cubic lattice. The system has the N^-xNyxN^ dimension and periodical boundary conditions are applied in all 
directions (Fig. y). We use N^ — Ny — 10 and Nz = 1000 in our simulation. The external electric field of the E 
strength is applied along the 2;-axis and the total sample width is L = dcNz. The site with intact cell is marked as 
occupied. We consider the idealised square pulse sequence with the pulse duration U, and the pulse repetition time 
At. 

The simulation procedure is done at each time step according to the following scheme: 

a) choose of the next occupied site within a lattice using the linear search procedure; 

b) determination of the number of occupied sites (n) among all of its six near- neighbours, determination of Pm and 
calculation of the transmembrane voltage m at a given cell with the help of Eqs. (|l|), (^ and (p^); 

c) calculation of the probability of a given cell destruction using Eq. (|l^) ; 

d) reiteration of the step (a) until all the occupied lattice sites appear to be tested. 

The total destruction degree P and values of ae and < /„(P, P„) > are calculated after each pulse. These time steps 
run until P reaches its asymptotic value. This model is a simple cellular automaton for the purposes of simulation of 
the heterogeneous material breakage kinetics. Indeed, at the step (b) we test the near-neighbour interior of a given 
cell and then we calculate the probability of destruction. 

The elevation of the temperature after each pulse was estimated with account to the Joule heating of material. The 
ohmic temperature increase AT{ti) after each pulse of ti duration was calculated using the following equation 

, , , dAT E'^CTeU , , 

^^(^^) = ^ = ^- (^0) 

We use in aU calculations: 4 = 10"^ m, d„, = 10"** m, ct^ = 0.1 - 1 S m^^, o-„, = 10"'' - 10"^ S mT^ , t^ = 10 - 100 
/iS and putting C ~ 3.93 kJ kg^-'^K"^ and p = 0.81xl0'^kg m""*, which values are characteristic for apples at 25°C 
(Losano,Urbicain & Rotstein,1979). 

III. MATERIALS AND METHODS 

A. Materials 

Freshly harvested apples of Golden Delicious variety were selected for investigation and stored at 4°C until required. 
In all the cases W was within 80-85%. We estimated the cell destruction degree as a ratio of effective conductivities 
measured before and after pulsed electric field treatment. The specific conductivities of samples, both initial (before 
treatment) and final (after the full treatment at highest voltages) were within ct^ = 0.004 — 0.008 S m~^ and Ca = 
0.1 — 0.2 S m~^, respectively. The specific conductivity of the apple juice extracted from the sample apples was within 



B. Methods 

Thin slices (6 ± 0.2 mm thickness and 45 ± 0.5 mm diameter) were cut from an apple pap. The conductivity was 
measured with LCR Meter HP 4284A (Hewlett Packard) for the thin apple slice samples and with Conductimetre 
HI8820N(Hanna Instruments, Portugal) for the apple juice samples at fr-equency 50Hz (this frequency was selected as 
an optimal in order to remove the influence of the polarising effects on electrodes and inside the samples) . Figure is 
a schematic representation of the experimental pulsed electric field treatment set-up. The temperature was recorded 
by the thermocouple THERMOCOAX type 2 (AB 25 NN, ±0.1°C). High voltage pulse generator, I500V-15A (Service 
Electronique UTC, France) allowed to vary ti within the interval of 10 — 1000 /iS (to precision ±2 /is), At within 
the interval of 1 — 100 ms (to precision ibO.lms) and N within the interval of 1 — 100000. All the experiments were 
repeated, at least, five times. Pulse protocols and all the output data (current, voltage, impedance, and temperature) 
were controlled with the data logger via Windows 95 software. 

IV. RESULTS AND DISCUSSION 

Figure |5| presents the experimental curves of relative conductivity (Ja/ <Jb versus number of pulses N for different 
values of the electric field strength E with ti = 100 /iS and Ai — 10ms. We have obtained practically similar results 
for (Ja/ub vs N dependencies at different values of At within the interval of 1 — lOOms. Hence, pulse repetition time 
does not influence our data and we have used Ai = 10ms in all the experiments. We have observed similar results for 
fa/cb versus Nti (equivalent time of electrical treatment) on ti variation. The kinetics of temperature evolution AT 
against N is presented in Fig. g[ 

We can conclude out of the data obtained that there exist, at least, two different stages of the material breakdown 
evolution. We have divided these stages conditionally by a horizontal dashed line in the Fig. |^. At the first stage, the 
sequential and the correlated breakdowns of cells develop in the system. The time of the first stage changes drastically 
with E increase. At low values of E we observe a very slow evolution of the material breakdown. The second stage of 
material breakdown (over the horizontal dashed line) flows more rapidly until terminated by an abrupt total dielectric 
breakdown of the material. At this moment we observe the overflow value of the out-of-limit current (15A) and stop 
the further treatment of material. We observe for each E a certain value of iVf,, which corresponds to the number 
of pulses at the moment of the total dielectric breakdown of material. Figure 0a presents N^ versus E (electric field 
strength) dependence. 

The insert in Fig. ra shows that practically in all the cases we observe approximately linear T versus N dependencies. 
Some deviation from this linear behaviour can be seen only for small E values {E — 0.2 kV cm~^) and at large 
N > W^; we can explain it by the heat exchange with the outside media. The slopes of the near-linear T versus N 
curves correspond to the mean values of AT(tj) = dAT/dN averaged over the total interval of N. They are shown 
in Fig. Rb for the different values of E. Here the solid line corresponds to the square-law AT{ti) versus E increase in 
accordance with relation 

AT{t,) = a + bE + cE'^, (21) 

which is consistent with Eq. (po|). Here a — 1.116x10^^, b = —1.118x10^^, c = 4.9463x10^^ are the values obtained 
from root mean square treatment of the experimental data for the interval E < 1 kV cm^^. At large values of field 
strength, E > 1 kV cm~^, we observed a considerable deviation from the parabolic law of Eq. ( po| ) in the dAT/dN 
versus E dependency. The reasons of such deviation are still unclear. 

The computer simulation results for the breakage degree (P) versus number of pulses (N) and effective conductivity 
(ae) versus number of pulses (iV) dependencies at the different values oi E, A = 0.1 and A = 0.01 are presented at Fig. 
m,b. We observe the obvious step- like behaviour of these P{N) and ae{N) curves describing the breakage kinetics. 
This behaviour reflects the "jamming" effects present in the systems under investigation and has a pure geometric 
origin. 

The "jamming" problem is very important and occurs in a number of situations, such as irreversible surface 
deposition of extended objects, random sequential adsorption (RSA problem), polymer physics problem or car-parking 
problem (J. W. Evans, 1993; Nielaba, Privman & Wang, 1993). The basic characteristic of this problem is the 
"jamming" coverage that depends on the type of the lattice and system dimensionality. In the case of deposition of 
the definite size particles the "jamming" limit is reached when it is impossible to place any further objects without 
overlapping the deposited before. The "jamming" effects in the dielectric breakdown problem have the following origin: 
when a cell gets punctured, there appear the bridging effects, which cause the abrupt fall of destruction probability 
for cells surrounding the punctured one. In the frame of our model, this bridging effect results from deviations of 



the local conductivity near the selected cell n (see, Eq. (JIJ)) against the average effective conductivity of the whole 
media (Te- 

In the RSA problem, the cell puncture is equivalent to the lattice site occupation, thus, occupation of definite 
site results in abrupt fall of the occupation probability for sites neighbouring the punctured one. The dielectric 
breakdown problem in such formulation is very similar to the RSA problem, except for successive "jamming" limits, 
which correspond to the site occupation with one, two, and et-cetera near-neighbours. We have carried out the Monte- 
Carlo simulation for the random sequential occupation of sites in the simple cubic lattice with j neighbour punctured 
cells in order to find the "jamming" limits P{j). We have obtained the following values for the lattice of 200x200x200 
dimension by averaging the calculation results over the 10 different initial configurations: P{0) — 0.305047 . . ., P(l) = 
0.42026 . . ., P(2) = 0.49368 . . ., P(3) = 0.544401 . . ., P(4) = 0.59321 . . ., P(5) = 0.64347 . . ., and P(6) = 1. 

The "jamming" effects in the dielectric breakdown problem result in existence of the saturation regimes in the 
breakdown kinetics. The first "jamming" step is clearly observed in the Fig. ^ for the small values of E/uo <1 at 
P — P(0) « 0.3. The next "jamming" steps are less pronounced and can be clearly observed only at small values 
of Au. The initial increase on the P versus N curves corresponds to the breakdown development without any near- 
neighbours bridging effects. At large values of E/uo > 1.5 we observe the jump on the next "jamming" steps (P(l), 
P(2)) with the subsequent saturation regime. 

We have observed the similar saturation regimes in experimental observations of the kinetics of dielectric breakage 
in thin apple slices (see Fig. ||), particularly, at small values oi E < 0.5 kV cm^^. The step-like behaviour is not 
so pronounced in experimental results as in computer simulation data presented in the Fig. [g; this difference results 
from the simulation model restrictions and, partially, from the lattice nature. 

Figure||c presents the curves of the calculated breakage degree P versus relative ohmic temperature increase Ar(°C) 
at different values of the electric strength E ior X = 0.1 and A = 0.01. These calculations were done using Eq. ([2^). 
These data allow us to understand the correlation existing between the achieved breakdown degree and the power 
consumption, which is proportional to AT. 

When the electric strength (E) values are low, the high destruction degree (P) can be achieved only on account of 
high power consumption and, correspondingly, with high overheat of the surrounding medium. However, the weakness 
of the electric field treatment is the possibility of a "jamming" regime, when the increase in power consumption fails 
to result in P increase, and poor control over the course of process. As far as it is difficult to choose precisely the 
required electric treatment mode for the given field strength values, here the transition to the mode of the overflowing 
out-of-limit current value occurs readily. In this case we stop the further treatment of a material and actually don't 
achieve the high values of P. 

Figure |^ presents the curves of the number of pulses in the total breakdown point (Nb) versus reduced electric 
strength [E/uo) for different values of Au/uo and A . We see that the character of Nb vs E curves depends on the 
values of Au/uo and A. The Nb increases with decrease of Au/uo at given E value. Unfortunately, we are unable 
to make more strict comparison between theoretical and experimental data, as far as we have no precise data for Uo, 
Am and A parameters present in computer model. However, the general character of calculated Nb vs E dependencies 
(Fig. 1^) correlates with experimental data presented in Fig. 0a. 

V. CONCLUSION 

The simplified dielectric breakage model based on effective media approximation is proposed. This approximation 
includes equations with the same morphology parameters as in percolation theory. The normalised transmembrane 
voltage / — u/dcE versus media breakage degree P dependence is obtained; this dependence is useful for estimating 
the transmembrane voltage and of the cell breakdown probability under the PEF treatment of cellular material. The 
most important feature of the proposed model is the existence of the "jamming" behaviour occurring in experimental 
observations of the biological tissue breakage. The different transitions corresponding to the different "jamming" steps 
are identified. The experimental study has yielded information about the material destruction degree and temperature 
elevation versus the time of PEF treatment at different values of field strength E in the interval of 0.2 — 2.2 kV cm"-'^. 
The proposed simulation model gives results consistent with the general trends observed in experimental breakage 
kinetics. The electric treatment is the most efficient at high values of the electric field strength E. However the 
efficiency of such electric treatment is limited by manifestations of the "jamming" effect, as well as by the processes of 
electric treatment upset. Reduction of E values allows controlling the treatment process more precisely and prevents 
from uncontrolled breakdown but it is accompanied with increase of electric power consumption. 
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FIG. 1. Normalised transmembrane voltage / — u/dcE versus media destruction degree P dependencies calculated for the 
MF approximation using Eqs. (H),(h) (curve 1), for the LC approximation using Eqs. (H), (112) (curve 2) and for the RSO model 
using results of simulation (curve 3). The dashed curve 4 corresponds to / dispersion A/ in RSO model. The calculations were 
performed at t = 2.3, Pp = 0.307, dc = lO"'^ m, dm = 10"* m, ac = 1 S m'^, am = W'^ S m'^, and A = 0.01. 



FIG. 2. Normalised transmembrane voltage u/dcE versus media destruction degree P at various n (the number of intact 
cells neighbouring a given cell) . The dashed line corresponds to the example of u/dcE behaviour with increase of the number 
of destroyed near-neighbours in the course of the media breakage. The calculations were performed at £ = 2.3, Pp = 0.307, 
dc = 10~^ m, dm = 10^* m, (Jc = 1 S m^^, am = 10^^ S m~^, and A = 0.01 for cells arranged in the sites of a simple cubic 
lattice. 



FIG. 3. The model of the cellular material structure used in the computer simulation. 

FIG. 4. Schematic representation of the experimental set-up used in the study of pulsed electric field treatment of the apple 
slices. 

FIG. 5. Relative conductivity (Ja/ob versus number of pulses A'^ at different values of the electric field strength E, ti = 100 
lis and At = 10 ms for thin apple slices at 25° C. 

FIG. 6. Relative ohmic temperature increase Ar(°C) versus number of pulses A'" at different values of the electric field 
strength E, ti — 100 fis and At = 10 ms for thin apple slices at 25° C. Insert shows the same data for small N values in linear 
co-ordinates T vs A*'. 

FIG. 7. Dependencies of the number of pulses in the total breakdown point Nb and of the mean ohmic temperature increase 
AT{ti) after each pulse of ti duration against the electric strength E a,t ti = 100 /xs and At = 10 ms for thin apple slices at 
25°C. 

FIG. 8. Calculated breakage degree P, effective conductivity ae versus number of pulses A'^ dependencies (a,b) and breakage 
degree P versus relative ohmic temperature increase AT (°C) dependencies (c) at different values of electric field strength E 
for A — 0.1 and A — 0.01. The calculation are performed at t = 2.3, Pp — 0.307, Au — 0.35, Uo = 1, ti — 100 ^s, dc — 10^"" m, 
dm = 10"** m, 0-c = 1 S m-^, am = 10"^ S m'^ ( A = 0.01), and a^ = 0.5 S m^S a„ = 5x10"^ S m^^ (A = 0.1). 

FIG. 9. Calculated number of pulses in total breakdown point Nb versus reduced electric strength E/uo for different values of 
Au/uo and A. The calculation are performed at t — 2.3, Pp — 0.307, Au = 0.35, Uo — 1, ti — 100 /is, dc = 10"^ m, dm ~ 10"* 
m,ac = lS m"\ am = 10"^ S m"^ ( A = 0.01), and ac = 0.5 S m"\ am = 5x10"^ S m"^ (A = 0.1). 



10 



This figure "fl.gif ' is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f8al.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f8bl.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f8cl.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f2.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f8a2.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f8b2.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f8c2.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f3.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f4.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f5.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f6.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f7.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



This figure "f9.gif" is available in "gif ' format from: 



http://arXiv.org/ps/cond-mat/9906 1 8 1 v 1 



